Statistics of Microstructure Formation in Martensitic Transitions Studied by a 
Random-Field Potts Model with Dipolar-like Interactions 
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We have developed a simple model for the study of a cubic to tetragonal martensitic transition, 
under athermal conditions, in systems with a certain amount of disorder. We have performed 
numerical simulations that allow for a statistical study of the dynamics of the transition when 
the system is driven from the high-temperature cubic phase to the low-temperature degenerate 
tetragonal phase. Our goal is to reveal the existence of kinetic constraints that arise from competition 
between the equivalent variants of the martensitic phase, and which prevent the system from reaching 
optimal final microstructures. 
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I. INTRODUCTION 



The kinetics of thernioelastic martensitic transitions 
(MT) is a challenging problem to which physicists and 
material scientists have devoted several decades of study. 
The complexity of the phenomenon arises from three 
main factors: (i) Firstly, the first-order character of the 
transition which, in many cases, implies the existence of 
high-energy barriers and associated long-lived metastable 
states. The word "athermal" has been used to refer to 
those MT where thermal fluctuations do not play any 
role, and thus only evolve through states when metasta- 
bility limits are reached. This, combined with the exis- 
tence of disorder in the systems, leads to history depen- 
dence and avalanche dynamics, (ii) The second impor- 
tant factor is the non-trivial relation between the symme- 
tries of the parent and product phases, especially in real 
3D systems. The product phase may appear in a num- 
ber of energetically equivalent variants which have the 
tendency to satisfy some " matching" conditions between 
each other and the parent phase, (iii) The third contribu- 
tion to the complexity of the phenomenon is the existence 
of long-range correlations due to the elastic nature of the 
problem. The effects of nucleation of a martensite do- 
main in a certain point of the system may influence the 
kinetics of the transition at a large distance from this 
point. The final microstructures observed in the marten- 
sitic phase are a consequence of the complex interplay of 
these three factors (i-iii) , and thus if one wants to have a 
full understanding of the phenomenon, it is not enough to 
reduce it to the problem of the optimization of a certain 
Hamiltonian. In order to understand this point of view, 
it is instructive to compare real microstructures observed 
over a large sample (see, for instance Fig.l in Ref. fT'| 
and Fig. 4 in Ref. with the perfect patterns calculated 
from elasticity theories (see, for instance Ref. [3|). The 
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definition of "kinetic constraint" in much simpler mod- 
els has been introduced [4, 5] to account for the fact that 
some low-energy states cannot be obtained due to the ab- 
sence of a possible pathway. Current theories are still far 
from being able to predict martensitic microstructures. 
Details such as sample size and shape, defects, annealing 
times, quenching rates, etc., are known to dramatically 
change the final state of the martensite (similarly to what 
happens for the domain structure of magnetic materials 
[li l3)- i&ct, it has been known for centuries that to 
obtain a "good" microstructure in a metallic alloy is an 
"art" reserved for the best smiths or metallurgists. 

Quantitative experimental information concerning the 
growth of the microstructure in martensitic transitions is 
difficult to be obtained [8] . This is because most observa- 
tion techniques, both direct imaging or scattering tech- 
niques (X-rays, neutrons and electrons microscopy), can- 
not be easily performed "in situ" during the transition. 
One must stop the transition in an intermediate state by 
quenching and perform a posterior sample treatment for 
the observation. Furthermore, in the case of direct imag- 
ing techniques it is difficult to extract numerical data 
from sequences of 2D micrographs of the system surface. 
There have been some "in situ" synchrotron radiation 
experiments that follow the growth of peaks associ- 
ated with martensitic structures, but the studies have 
mostly concentrated on understanding the precursor ef- 
fects before the real transition starts. Some interesting 
phenomena have been reported during the transition. In 
some cases the formation of intermediate "phases" with 
only short-range order and with reciprocal space vectors, 
different from the final position of the martensitic peak 
structure have been observed 0]. In addition, a recent 
work has quantified the fact that microstructures exhibit 
imperfect self-accommodation due to the coalescence of 
the variants |]|. 

Among the theoretical efforts to understand mi- 
crostructure formation, one should mention many contin- 
uum models, derived from elasticity theory [H, IH, 
[isl . [T^ . [Tsl . [l6j . Despite advances in computing power, 
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the complexity of the models (especially the long-range 
character of the interactions) is such that it is difficult to 
perform a large number of simulations and a statistical 
analysis of microstructures. In many cases the models 
have only been applied to 2D unrealistic crystals, while 
in others the number of coexisting variants has been re- 
duced artificially or by considering a symmetry breaking 
external stress. 

In this work we will focus on the cubic to tetragonal 
MT. This kind of transition occurs, for instance, in a 
number of binary metallic alloys for different ranges of 
concentrations. Among others Q we recall InTl, InPb, 
NiAl, FePt and FePd ^9J. Different microstructures have 
been observed in these transitions. A common basic fea- 
ture is the clear tendency for the variants to grow forming 
twins. These are domains of different variants that share 
a common interface oriented in a certain preferred direc- 
tion. By this mechanism, elastic distortions of the lattice 
are minimized. For the cubic to tetragonal transition of 
interest here, the twin boundaries are planes perpendic- 
ular to the (110) direction (or any of the 6 equivalent 
ones) of the cubic phase. The overall microstructures, 
nevertheless, are a combination of the different twins and 
usually look quite random. 

In a previous paper we introduced a Random Field 
Potts Model with truncated dipolar interaction for 
"athermal" first-order phase transitions between phases 
with any change of symmetry [l3| . The model was based 
on the T = RFIM with metastable dynamics, which 
was originally proposed for the study of field induced 
athermal transitions between two ferromagnetic symmet- 
rically equivalent phases. The proposed modifications 
allowed, for instance, the following to be studied: transi- 
tions from a non-degenerate phase to a phase with three 
different variants, the influence of the amount of disorder 
on the dynamics and how the parameters of the dipolar 
interaction term affect the final microstructure. 

In the present study we are going to go one step for- 
ward. We will adapt the model to the symmetries of 
a cubic to tetragonal MT. We will extend the range of 
the dipolar interaction in order to obtain the twin mi- 
crostructures observed in such a MT transition. We have 
found that it is enough to include up to 4th nearest neigh- 
bors. By truncating long-range interactions we acceler- 
ate our simulations, so that we can perform a systematic 
study of the Hamiltonian parameters and statistics over 
disorder configurations. 

The paper is organized as follows: in Section |TT] we in- 
troduce the Hamiltonian of our model, and discuss the 
truncation of the dipolar term. We also give details 
about the dynamics and the simulations. In Sec lIII A] 
we present a ground state analysis of the martensitic mi- 
crostructures that optimize the energy. This will help us 
in choosing the interesting range of Hamiltonian param- 
eters. In Sec. nil B I we present some of the microstruc- 
tures obtained by dynamically driving the system from 
the parent to the fully saturated martensitic phase. The 
results are compared with the previous ground-state con- 



figurations. In Sec lIII C] we discuss the shape of the hys- 
teresis cycles and the effect of varying the Hamiltonian 
parameters (disorder and dipolar intensity) and the sam- 
ple size. A more quantitative analysis of the microstruc- 
tures is presented in Sec lIII Dl and IIIIEI were we show 
the system evolution along the hysteresis cycle in real 
and Fourier space, respectively. Finally in Sec. IIVI we 
both summarize and conclude. 



II. MODEL 

Let us consider a simple cubic lattice of size N — 
L X L X L, with lattice parameter a and periodic bound- 
ary conditions. At each lattice site we define a variable 
which can take four different values that we will call 0, 
i, y and z. As in our previous work [l3|, we have cho- 
sen to represent our variables by considering a vector Si 
{i = l,...iV), having three components: we will indi- 
cate the four possible values as = (0, 0, 0), x = (1, 0, 0), 
y = (0,1,0) and z = (0,0,1). These variables can be 
interpreted as an elementary domain in the cubic austen- 
ite phase (vector 0), or in the three possible variants of 
a tetragonal martensitic phase {x, y and z). The in- 
teraction between the variants can be described by the 
following Hamiltonian: 
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FIG. 1: (Color on line) Schematic representation of the near- 
est neighbors in a cubic lattice, indicating the NN (continuous 
line) and 4N (dashed line) positions. 

The first contribution in the Hamiltonian is a Potts 
exchange term between nearest neighbors (NN) which 
accounts for the bulk energy favoring the growth of sin- 
gle phase domains. The second is the long-range dipolar 
interaction. The absolute value taken in the dipolar term 
guarantees the rotational invariance that is expected in 
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the systems that we want to model. The third Hamilto- 
nian term represents the interaction between the external 
driving field H and the order parameter M — (-^i)^- 
Note that the vector 0, corresponding to an elementary 
domain in the austenite phase, has absolute value zero 
and thus does not contribute to M . On the other hand, 
all the three variants of the martensitic phase are here 
represented by vectors with modulus one. Consequently, 
M/N can be read as the normalized transformed frac- 
tion of the sample. Thus, the field H drives the system 
from the cubic phase m = (when H — —oo) to a multi- 
variant tetragonal phase with m ~ 1 (when H = oo). It 
therefore plays the role of the temperature for such an 
"athermal" phase transition. 

The last term of the Hamiltonian accounts for the in- 
teraction with the quenched disorder (impurities, dislo- 
cations, vacancies, etc.), which is always present in real 
materials. Its precise form will be discussed later on in 
this section. 

As explained in the introduction, a first version of this 
model considered only nearest-neighbor dipolar interac- 
tions [l3|. This produced anisotropic domains of the 
different variants, but not the correct microstructures 
observed in MT. In the present work we will introduce 
higher order dipolar interactions. As will be seen, it is 
enough to include the interaction between fourth near- 
est neighbors (4N), at a distance 2a from the reference 
elementary domain (see Fig. [T]). It will be the interplay 
between the interactions to NN and to 4N (which are 
placed along the same spatial direction) that generates 
the convenient microstructures^n a way similar to what 
happens in the ANNNI model |i3 • 

Although we have also studied the second and third 
nearest-neighbors interactions, we have found out that 
these terms are not crucial for the microstructurc for- 
mation and do not add any new physics. Therefore, for 
simplicity in most of the paper they have been neglected, 
except where indicated otherwise. Moreover, we have al- 
lowed the coefficients A/ |ry|'^ that control the decay of 
the dipolar interaction to have more freedom. Since we 
are truncating and including only two terms, there is no 
need to keep the coefficients as in the original model. We 
have considered two general parameters Ai and A2 multi- 
plying the dipolar interaction to NN and to 4N nearest- 
neighbors. Thus, after these modification, our Hamilto- 
nian reads: 

NN NN 

n^-J2 ^is..' ^i) + Ai 5^ Ks'. • h/m ■ %)l 

<ij> ij 
4N N 

+ A2 5^ 1(5. • u,){s, ■ h,)\ -HY,is^r + nd^s, (2) 

ij i 

where the second sum runs over NN and the third only 
over the fourth nearest-neighbors. The interaction with 
quenched disorder has also been simplified from what was 
originally proposed in Ref. Here, in fact, we will not 
focus on the study of the influence of disorder, but we 



need a certain amount of it in order to favor nucleation 
of the martensitic phase in all the equivalent variants. 
Thus we have considered: 

N 

Ti-dis = — o" ^ 9i- Si, (3) 

i 

where the quenched random fields cji are independent and 
identically distributed Gaussian random variables, with 
zero mean and unitary variance. The parameter a models 
the strength of the disorder. Note that after the simplifi- 
cations, the Hamiltonian parameters are only three: Ai, 
A2, and (T. 

For the study of athermal behavior, apart from the 
Hamiltonian, one must choose a certain dynamics. We 
chose to implement the extremal update dynamics p7| : 
starting from a saturated system configuration, corre- 
sponding to m = at = — 00, we increase (decrease) 
the field by small steps AT?. At each field step we check 
the contribution to the energy of each Si. When a vari- 
able is found to decrease the total energy by varying to a 
new state, we change it to the value that gives the most 
negative energy change. In this way we guarantee that 
the system will reach the same state for the same applied 
field values, independently of the AH value. The states 
reached are equivalent both macroscopically (same trans- 
formed fraction m) and microscopically (same values of 
Si for all the i values). The simulations presented below 
were obtained by choosing AH = 0.05. 

III. RESULTS 
A. Ground state microstructures 

Starting from the experimental observations from cu- 
bic to tetragonal transitions, our goal would be the 
growth of the so-called twin variants. These correspond 
to regions of alternating x and y (e.g.) domains sepa- 
rated by interfaces parallel to the z axis and forming an 
angle of 45° with the x and y axis. Equivalent structures 
will be the x-z twins separated by interfaces parallel to 
the y axis and forming 45° with the x and z axis, and 
with the y-z twins separated by interfaces parallel to the 
z axis and forming 45° with the x and y axis. 

In addition, we would be interested in controlling the 
width w of such twins. If one truncates the dipolar term 
to only first and second nearest-neighbors, it is not pos- 
sible to stabilize such twin variants, except for the case 
with width w = 1, which in fact corresponds to do- 
mains displaying a chessboard structure as illustrated in 
Fig. [2l^a). Due to lattice geometry, third-nearest neigh- 
bor interaction can not lead to wider structures and it is 
only by including 4N terms that we can obtain a. w — 2 
structure, which is represented in Fig. Hfb). 

In order to study the range of parameters in which 
the different microstructures occur, we have performed a 
comparison of the energies of the different configurations 
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shown in Fig. [2l with w = 1,2,3, 6, cxd. The analyti- 
cal results presented here are general for any system size 
provided that the proposed configurations are compatible 
with periodic boundary conditions. 



From such energy functions one can obtain the phase 
diagram that indicates the boundaries between the struc- 
tures that minimize the energy. This is shown in FiglS] 
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FIG. 2: (Color online) System configurations of y ~ z twins 
with widths w = 1,2,3,6,00, plotted for a system with size 
L = 12. 

The values of the energy as a function of the model 
parameters (in the absence of disorder) of each of the 
configurations considered above are given by: 
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FIG. 3: (Color online) Phase diagram in the A1-A2 plane. 
The solid lines represent the interfaces between the phases 
obtained from Eqs|4l 

As can be seen, the Eqs. |3] allow the existence of 
ground-state twinned structures w = 1, w — 2, and 
w — 00, while no other phases can exist. Although the 
w = 2 structure is good enough for our objective of study- 
ing kinetic constraints, we have also analyzed the mod- 
ifications one should perform on the model in order to 
obtain ground-state twin structures with larger widths. 
The guess is that higher order dipolar terms would lead 
to the stabilization of the wider w domains. To confirm 
this guess, let us once again consider a. d ~ 12 system, 
and add the seventh nearest-neighbor term (that is, the 
neighbors at a distance of 3a from the reference elemen- 
tary domain). The energies of the various phases are 
now: 

Eoo = -3N + NXi+N\2+NX3 
El = -N + NX2 

E2 = ^2N + N^Xi + ^X3 

E3 - -IiV + iV^Ai+iviA2. (5) 



El = -N + NX2 
E2 = -2N + N^Xi 

E3 = -'^N + N^Xi + Nh2 

47 97 

-Boo = -3N + NX1+NX2. (4) 



An example of a projection of the phase diagram is shown 
in Fig. m where it can be seen that the w — 3 phase is 
now stable. The borders between the phases are calcu- 
lated with Eqs. [H 

For the remainder of the paper we will concentrate on 
the A3 = case, and in the region where w = 2. Most 
of the results presented will correspond to Ai — 10 and 
A2 = 20. 
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FIG. 4; (Color online) Phase diagram in the plane A1-A3 for 
A2 = 2. 



B. Kinetically obtained microstructures 



In Fig (Ha) we show an example of microstructures ob- 
tained by the numerical simulations of the Hamiltonian 
(Eq. [2|) using athermal dynamics. We represent only 
three perpendicular projections of our 3D calculation. 
This example thus only supplies qualitative information. 
In section llll El we will perform a Fourier analysis in order 
to gain some global information. 

From the Figure (Hla) it is clear that some local struc- 
tures of the kind w — 2 arise, but that we are very far 
from the ground-states computed previously (see Figl2]). 
The global structure is more complex and displays mixing 
of the three possible twin domains. Taking into account 
that there are 3 twin pairs and that each of them can oc- 
cur in 2 different equivalent orientations and that, given 
our underlying discrete lattice, there are 4 possible trans- 
lations (for our w = 2 case) , the number of different twin 
domains is 24. 

Another interesting example is given is Fig. [5fb). This 
figure corresponds to the saturation configuration ob- 
tained by the application of a very high external com- 
pression in the z direction, disfavoring the associated z 
variant. In this case, the formation of diagonal w = 2 do- 
mains in the x — y plane is already very clear at first sight, 
confirming the importance of kinetic constraints arising 
due to the coexistence of the three equivalent variants. 
It is interesting too to analyze the microstructure of the 
X — z and y — z planes. In this case, quite elongated 
domains formed by vertical alternating stripes with pe- 
riodicity 2 arise, as indicated in the figure. This kind 
of domains have been observed experimentally in a iron- 
palladium ferromagnetic shape-memory alloy, displaying 
a cubic to tetragonal martensitic transition f^. 
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FIG. 5: (Color online) Examples of two microstructures at 
saturation for Ai = 10, A2 = 20, = 2. and L = 36. In 
(a), the three variant saturation microstructure, in (b), the 
stressed configuration obtained by the application of a high 
external compression in the z direction. White lines in the 
yz planes are guides to the eye for the identification of the 
regions of twinned domains. 



C. Hysteresis 

In this section we study the effect of the system size 
and of the Hamiltonian parameters on the hysteresis cy- 
cles. As will be seen in all of the figures of this section, the 
two branches of the loops are strongly asymmetric, due 
to the physical difference between the process of transfor- 
mation from a one- variant phase to a three- variant phase 
(ascending branch) and vice-versa in the descending case. 

In Fig. [HI we show some examples of hysteresis loops 
obtained by varying the amount of disorder a. The loops 
have a tendency to display plateaus, which are smooth 
when disorder increases. An additional effect of disorder 
is also to tilt the hysteresis loops. 

In Figs. [21 we present examples of the effect of the 
dipolar parameters Ai and A2. Within the region of sta- 
bility of the w = 2 phase, the plateaus are always present, 
provided that disorder is not too high. 
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FIG. 6: (Color online) Hysteresis cycles for a — 0.01; 2; 5, for 
system parameters L = 24, Ai — 10, and A2 = 20. 



FIG. 8: (Color online) Hysteresis loops for L = 12, 24, 36, 48, 
for system parameters Ai — 10, A2 = 20 and a = 2. 
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FIG. 7: (Color online) Hysteresis cycles for (Ai,A2) = 
(1, 2); (5, 10); (10, 20), for system parameters L = 24, and 
(7 = 2 (a) and a = 0.5 (b). 



Finally, in Fig. [8] we study the effect of system size on 
the hysteresis loops. The loops do not show a strong size 
dependence, and the plateaus are clearly not a finite-size 
effect. In the following we will restrict our simulations to 
the affordable intermediate size L = 36. 

All the hysteresis cycles that we have shown in this sec- 
tion display qualitatively the same behavior, consisting 
of a certain number (2-3) of plateaus, where the trans- 
formed fraction remains constant in a range of H values. 
This feature suggests that the transformation of the 
phase has a tendency to split into steps, the first involv- 
ing the transformation of more than 50% of the lattice. 
In the next sections we will try to elucidate the physical 
reasons behind this tendency. 



D. Real space analysis of the transition dynamics 

In order to investigate the nature of the intermediate 
steps of the hysteresis cycles we can analyze the system 
configurations in real space at various points of the loops, 
as indicated in Fig[9^a). We also show the derivative 
dm/dH which better reflects the activity of the transition 
and reveals the importance of the first transformation 
step compared to the others. These kinds of activity 
curves can be compared with the experimental data from 
calorimetry or acoustic emission activity. 




FIG. 9: (Color online) Hysteresis ascendant loop branch (a) 
and (b)field derivative of the normalized transformed frac- 
tion dm/dH as a function of the field, for L — 36, Ai = 10, 
A2 = 20, and cr = 2. The points marked with numbers 
1, 2, 3, 4, 5, 6, 7, 8 correspond to snapshots in Fig. 1 101 and have 
m = 0.125; 0.25; 0.375; 0.5; 0.625; 0.75; 0.875; 1 respectively. 

The microstructures on each indicated point of the as- 
cending branch are shown in Fig. 1101 A first inspection 
reveals that the transition basically proceeds by nucle- 
ation events. There is some growth, but almost no coars- 
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FIG. 10: (Color online) Snapshots of real space configurations corresponding to the field loop points indicated in Fig. [5] 



ening since the possibilities for two domains to match are 
small. One can also see that re-transformation events 
are common. The observed step in the hysteresis loop 
is clearly associated with the fact that martensite do- 
mains percolate in the sample. At this point the system 
needs an extra increase of the driving force to proceed. 
The subsequent domains will grow surrounded by a com- 
plex boundary of the martensitic phase. Therefore, the 
growth dynamics will be very different from the nucle- 
ation of the first domains in a cubic phase environment. 
In the second part, there are many re-transformation 
events until the final saturated microstructure is ob- 
tained. After a more careful study of the snapshots one 
can infer that in the first step of the transformation, 
corresponding to the snapshots from 1 to 6, the system 
chooses to develop oblate (disk-like) domains. For in- 
stance, one can see in snapshot number 3 that the i-phase 
(black) tends to be distributed in extended regions in the 
y — z plane, whereas it shows elongated domains perpen- 
dicular to the X-axis on the x — y and x — z planes. Such 
structures correspond to the structures that optimize the 
NN dipolar interaction term (Ai), as we studied in a pre- 
liminary work p7j . In a second step, when the external 
field is comparable to the value of the parameter A2 , the 
system is able to develop the w = 2 structures with the 
desired orientation, and one starts to see regions of twin 
domains with the 45° orientation on the different planes. 
Of course, they are very much influenced by the domains 
that have already grown. This interpretation can even 



be argued by confronting the curves for different values 
of Ai and A2 in Fig. [71 For a better quantitative un- 
derstanding of this effect, it is convenient to perform the 
analysis of the microstructures in Fourier space. 



E. Fourier analysis 

A quite natural way to analyze the system microstruc- 
ture is to study the Fourier Transform (FT) of the density 
function. In fact, as it is well known, the square modu- 
lus of this quantity represents the intensity that can be 
obtained from scattering experiments. In our case, at 
saturation we have the coexistence of three variants, x, 
y, and z, and thus it is useful to separately consider three 
density matrices Px{x, y, z), Py{x, y, z), and Pz{x, y, z). In 
analogy we can define the density matrix of the phase 
PQ{x,y,z), and transform the four matrices individually. 
Thus each variant a [a — x, y, z) is associated with a 
scattered intensity: 



E 

b.c.d=l 



Pa{b, c, d)e' ^ 



(K^b+KyC+K^d) 



where r — xi -\- yj -\- zk — bai -\- caj + dak is the generic 
position vector, and K^, Ky, = 0, L — 1. The elements 
of these four density matrices can take the values 1 or 0, 
representing, respectively, the occupation or the vacancy 
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FIG. 11: Positions of the peaks of the scattered intensity 
for the ideal case with w = 2, for a system with L — 
36 and without disorder. The configuration is obtained as a 
superposition of the six possible orientations of the w = 2 
stripes. 



of a lattice site by a variant of a given color. The presence 
of disorder can locally favors one or the other variant, 
but since the disorder is randomly extracted, we expect 
variants behavior to be equivalent on average. We thus 
restrict ourselves to studying just one variant, e.g. the x 
one, without loss of generality. In the simulations we have 
verified that the color equivalence is effectively fulfilled. 
In the following we will address only this case. 

In Fig[Tl] we show the position of the peaks of the 
squared Fourier Transform of the density of the x vari- 
ant in the ideal case w = 2. We observe that Bragg peaks 
occur at the position K — (L/4, L/4, 0) and in symmetric 
equivalent positions. In fact, if one performs the Fourier 
transform of the density of i-phase corresponding, for 
instance, to the snapshot in Fig[21 one obtains only two 
of the eight Bragg reflections shown in FiglTTl Neverthe- 
less, for symmetry reasons one has to consider the two 
possible 45° orientations and also that the phase x would 
be present not only in the x — y twin, but also in the x — z 
twin. Translations of the configuration with the same do- 
main orientation give the same contribution. Therefore, 
we shall expect the 8 Bragg reflections shown in FigfTTl 
with equal intensity 4(^/4, L/4, 0) = l/6{N^/8), where 
the contribution 7V^/8 is given by the intensity of the two 
secondary peaks that arise for each of the four ideal con- 
figurations where the variant x is present, and the factor 
1/6 is due to the average over all the six possible ideal 
configurations. The peak at the origin stands for the 
fraction of the system which is covered by the x phase, 
and thus its intensity is given by /^(0, 0, 0) = 4/6(iV^/4), 
where again the factor 4/6 is due to the average over the 



six possibilities. In a real experimental case, the widths of 
the twin domains will be much bigger than the w = 2 lat- 
tice spacing and, therefore the growth of the martensitic 
tetragonal domains will be revealed by satellite peaks 
occurring close to the cubic Bragg reflections in the di- 
rections (1,1,0) (and equivalent ones). 

In Fig. [T^l we show the Fourier analysis of the mi- 
crostructures obtained after saturating the system by the 
numerical simulations. The intensity is averaged over 100 
realizations of disorder. We see that the peaks appear at 
the positions expected from the ground-state structures, 
but the widths of the peaks are quite large. 
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FIG. 12: Scattered intensity (linear scale): Fourier trans- 
form of the X particle density for a sample with L — 36, a = 2 
Ai — 10, and A2 = 20. Data are averaged over 100 realizations 
of the disorder. The K = (0, 0, 0) peak has been suppressed 
in order to enhance the contrast. 

In analogy with the study performed in real space, we 
can study the evolution of the Fourier transform along 
the hysteresis cycles. This is presented in Fig. [T31 As we 
have already noticed, the system evolution begins with 
a tendency to form and grow disk-like domains which 
correspond to peaks &t K = (L/4, 0, 0) (and equivalent). 
After the first plateau the intensity is transferred to the 
final positions ai K = {L / A, L / A, 0) (and equivalent). 

In order to gain a quantitative insight into the dynam- 
ics of the transition, we have analyzed the evolution of 
the scattered intensity corresponding to some significant 
vectors of the Fourier space. In Fig. [HT aV we present 
the evolution of the diagonal — Ky with = 0. As 
can be seen, the peak is slightly shifted from the ideal po- 
sition (L/4, L/4,0) (and equivalent). In particular, the 
position is closer to the origin, meaning that the forma- 
tion of some thinner structure is affecting the average. 
This effect is due to the presence of the disorder and 
of the dynamic constraints. In Fig lMf b). we present the 
evolution of the intensity as a function of the transformed 
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FIG. 13: Scattered intensity Is:{K) (linear scale) for g = 2, \i = 10, A2 = 20, and L = 36. Data are averaged over 100 disorder 
realizations. The snapshots correspond to the points indicated in Fig. [9] 1) m = 0.125, 2) m = 0.25, 3) m = 0.375, 4) m = 0.5, 
5) m = 0.625, 6) m = 0.75, 7) m = 0.875, and 8) m = 1. The K = (0, 0, 0) peaks have been suppressed in order to enhance the 
contrast. 
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fraction m. We choose m and not H as the independent 
variable since, as can be seen in Figs. [7] and El the 
jumps in the hysteresis cycles are very sharp. As can be 
seen, the K = (0,0,0) and the K = (L/4,L/4,0) val- 
ues grow monotonically with m as one could expect, and 
the intensity at K ^ (L/4, L/A, 0) and K = {L/2, L/2, 0) 
tends to disappear at saturation. The scattered inten- 
sity in this latter position is a measure of the disorder 
in the system, since in the ideal configuration the corre- 
sponding value of Ix is zero. In the figure we confront 
the K = (0,0,0) and the K = (L/4,L/4,0) values with 
the related theoretical values at saturation. The value of 
the K = (0, 0, 0) peak for m = 1 is in perfect agreement 
with the expected one, while the K — (L/4, L/4, 0), even 
if increasing, as it should do in the presence of w — 2 
structures, keeps far from the ideal value. Furthermore, 
the predominance of the w = 2 structures arises only in 
the last part of the hysteresis branch, where the reduced 
transformed fraction goes from m = 0.875 to m = 1. 
These effects are due to the presence of disorder and of 
the kinetic constraints associated with the presence of 
many competing growing domains. 
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IV. CONCLUSIONS 



We have presented a model for the statistical study 
of microstructures in martensitic phase transitions. We 
have focused on the cubic to tetragonal case. The model 
is not aimed at exactly reproducing the microstructure 
details but only to provide a statistical picture of the dif- 
ferent phenomena that may occur during the transition 
and that will constrain the final martensitic state. 

Our model allows some dynamical aspects of the MT 
that occur due to the existence of kinetic constraints to 
be studied. Among others, we have shown the existence 
of domain re-transformation, the fact that the transi- 
tion may proceed through several steps, and the exis- 
tence of intermediate structures that minimize only the 
short-range term of the dipolar interaction. The Fourier 
transform of the obtained microstructures has enabled 
a quantitative analysis of such phenomena after taking 
averages over disorder configurations. Although some of 
the observed phenomena may not be totally realistic due 
to the cut-off that we have introduced in the long-range 
interaction, we have shown that the model is good enough 
to exhibit kinetic constraints. The possibility for extend- 
ing the interactions to real long-range interactions is still 
open and, although it will require much longer computing 
times, should not be excluded for a future work. 

In our opinion, the presented results quite clearly 
justify the importance of a statistical approach to the 
martensitic phase transition. We finally propose that, be- 
sides the real imaging experiments, many more "in situ" 
scattering experiments will be desirable in order to gain 
understanding into the dynamics of microstructure for- 
mation in martensites. 
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FIG. 14: (a) Intensity h for m = 

0.125; 0.25; 0.375; 0.5; 0.675; 0.75; 0.875; 1 for the ascen- 
dant branch on the Kx = Ky diagonal, with = 0. (b) 
Values of the scattered intensity as a function of the reduced 
transformed fraction, for K = (0,0,0), K = (L/2, L/2, 0), 
K = (L/4, L/4, 0), and K = (L/4, 0,0). Dot-dashed and 
dashed line represent the theoretical values of satura- 
tion for ideal configurations of the K — (0, 0, 0) and the 
K = (L/4, L/4, 0) peaks. The values are averaged over 100 
disorder configurations, for system parameters cr = 2, p = 0, 
Ai = 10, A2 = 20, and L = 36. 
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